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We derive a supersymmetric renormalization group (RG) equation for the scale-dependent super- 
potential of the supersymmetric 0{N) model in three dimensions. For a supersymmetric optimized 
regulator function we solve the RG equation for the superpotential exactly in the large-A'' limit. 
The fixed-point solutions are classified by an exactly marginal coupling. In the weakly coupled 
regime there exists a unique fixed point solution, for intermediate couplings we find two separate 
fixed point solutions and in the strong coupling regime no globally defined fixed-point potentials 
exist. We determine the exact critical exponents both for the superpotential and the associated 
scalar potential. Finally we relate the high-temperature limit of the four-dimensional theory to the 
Wilson-Fisher fixed point of the purely scalar theory. 
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I. INTRODUCTION 

Fixed points of the renormalization group (RG) play 
a fundamental role in statistical physics and quantum 
field theory [U [2]. Infrared (IR) fixed points dominate 
the long-distance behavior of correlation functions and 
are relevant for the understanding of continuous phase 
transitions and universal scaling laws ^. Ultraviolet 
(UV) fixed points control the short-distance behavior of 
quantum field theories. It is vifidely believed that the 
existence of an UV fixed point is mandatory for a def- 
inition of quantum field theory on a microscopic level, 
e. g. asymptotic freedom of QCD or asymptotic safety of 
gravity [H [S] . In general, the fixed point structure of a 
given theory depends on its field content, the spacetime 
dimensionality, the long-range or short-range nature of 
its interactions and the symmetries of the action. 

Scalar field theories with a global 0{N) symmetry pro- 
vide an important testing ground for fixed point stud- 
ies. In three dimensions the (0^)^ theory displays a non- 
trivial IR fixed point which determines the second-order 
phase transition between an 0{N) symmetric and the 
symmetry broken phase as realized in many physical sys- 
tems ranging from entangled polymers and water to fer- 
romagnets or QCD with two massless flavors of quarks 
[3i_6^. The (0^)'^ theory also displays a line of flrst-order 
phase transitions whose end point, in the limit of many 
scalar fields, quahfies as an UV fixed point [71 [5]. 

Supersymmetry represents the global symmetry which 
relates bosonic to fcrmionic degrees of freedom. Su- 
persymmetric theories are important candidates for ex- 
tensions of the Standard Model. It is important to 
understand how the fixed-point structure of a non- 
supersymmetric theory differs from that of its supersym- 
metric extension, both in view of the IR and the UV 
behavior of the theory. 

In this paper, we study fixed points of supersymmetric 
0{N) models which consist of an A^-component scalar 
field coupled to N Majorana fermions. We employ non- 
perturbative renormalization group methods a la Wilson, 



based on the integrating-out of momentum modes from 
a path-integral representation of the theory [51-ITT]. A 
particular strength of this continuum method is its flex- 
iblity, allowing for the study of theories with strong cor- 
relations and large couplings. Furthermore, optimization 
techniques are available to control the physics content 
within systematic approximations |12H14j . In the past, 
this method has been successfully employed for the study 
of critical phenomena in a variety of settings includ- 
ing scalar theories, fermions, gauge theories and gravity 
[T5HM| . It has recently been extended to include su- 
persymmetric theories p5ti36j . Our prime interest here 
concerns the limit of many scalar fields 1/N — >■ 0, where 
effects induced by the fields' anomalous dimensions are 
suppressed and a local potential approximation (LPA) 
becomes exact. Then full analytical fixed point results 
are obtained for the fixed points in the supersymmetric 
theory, allowing for a complete analytical understanding 
of the theory, analogous to the purely scalar theory [371 - 

ES]. 

Supersymmetric 0{N) models have previously been in- 
vestigated with Dyson-Schwinger equations [40j and with 
the large-A^ expansion [HI 22] • The three-dimensional 
theory has also been studied at finite temperatures, 
where supersymmetry is softly broken [331 |33]- The 
model has a peculiar phase structure concerning the 
breaking of the 0{N) symmetry: Additionally to the 
normal phases with a broken and an unbroken symme- 
try a phase with two 0{N) symmetric ground states 
and a phase with one symmetric and one non-symmetric 
ground state have been found. In addition, there exists a 
supersymmetric analogue of the Bardeen-Moshe-Bander 
phenomenon [71 . The fate of this phenomenon at finite 
A^ remains yet to be resolved [151H7] . 

The paper is organized as follows: First we introduce 
the supersymmetric 0{N) model (Sec. TO and derive 
the non-perturbat ive flow equation for the superpoten- 
tial in LPA (Sec. |III[). We then solve this equation an- 



alytically in the large-A^ limit (Sec. IIIC) and analyze 
the resulting fixed-point solutions (Sec. IV). We com- 



pute the universal scaling exponents and compare our re- 



suits with those in the non-supersymmetric theory with- 
out fermions (Sec. Ivl). We close with a discussion of our 
results (Sec. |VI[ ). Our conventions and a derivation of 
supersymmetric flow equations in superspace is found in 
the appendix. 



II. SUPERSYMMETRY 

In this section, we recall the definition of three- 
dimensional supersymmetric 0{N) models, which are 
built from N real superfields 

<^'{x,0)^(l)' + e^//{x) + heF\x), i = i,...,N. (1) 

Each component of the superfield contains a real scalar 
field, a two-component Majorana spinor field and a real 
auxiliary field, $' ~ ((/>%'0% i^*). We shall use a Majo- 
rana representation with imaginary 7-matrices {7''} — 
{(T2,icr3,Z(Ti}. Then the metric in {7'^, 7"^} = 27^'"^ takes 
the form rj^^ = diag(l, — 1, — 1). A Majorana spinor is 
real in this representation and V' = {iip2, —itpi)- The 
supersymmetry variation of the superfield is generated 
by the supercharge Q via ^e$' — zeQ$% where the ex- 
plicit form of the supercharge and further conventions 
are collected in appendix [A] To construct a supersym- 
metric invariant action we note that the i^-term in the 
expansion (fTj) transforms under supersymmetry transfor- 
mations into a spacetime divergence such that its space- 
time integral is invariant. 

In order to define an 0{N) symmetric, supersymmetric 
action we introduce the supercovariant derivatives 

^+0e and V^----ie0, 
06 aO 



V 



(2) 



which anticommute with the supercharges and thus map 
superfields into superfields. Since the theory should be 
0{N) invariant, the superpotential only depends on the 
invariant composite superfield i? = ^^^^i. In component 



form, it reads 



R=g+{04'^W 



'F, - -^V, 



(3) 



where the quantity g = ^'f'^'t'i has been introduced. The 
starting point for further investigations will be the su- 
persymmetric action 



S 



1 



-<i>'VV<!>, + 2NW 



(4) 



which contains a kinetic term with supercovariant Lapla- 
cian VV as well as an interaction term, given by the su- 
perpotential W. We have already rescaled the fields and 
the superpotential with N. An expansion in component 
fields yields the Lagrangian density 



vr (I) *'*-»'" (I) 



'(I)*- 



w ^^ 



F' 



2N 



(5) 



where primes denote derivatives with respect to g/N. 
Eliminating the auxiliary field F^ by its algebraic equa- 
tion of motion, F* = ~W'{g/N)(j)'', yields the on-shell 
Lagrangian density 



I i -■ 

-Con = --0''n0, - -ip'^ 



-^^"'11) 



-w" 



N 



2N 



(6) 



From ^ we conclude that the potential for the bosonic 
field follows from the superpotential W via 



Vi-g) = -gW-{l) 



(7) 



Note that for a polynomial superpotential W{g/N) which 
for large g tends to W '^ g" we do not expect supersym- 
metry breaking in our non-perturbative renormalization 
group studies. 



III. RENORMALIZATION GROUP 

A. Supersymmetric flow^s 

In order to analyze the phase transition and the low- 
energy behavior of supersymmetric sigma models we 
resort to Wilsonian renormalization group techniques. 
Specifically, we adopt the framework of the effective av- 
erage action based on the infinitesimal integrating-out of 
degrees of freedom with momenta q"^ larger than some 
infrared momentum scale k^. In consequence, the effec- 
tive action becomes a scale-dependent effective action Tk 
which interpolates between the microscopic action S in 
the UV and the full quantum effective action in the IR, 
where fc — > 0. The scale-dependence of F^ is given by an 
exact functional differential equation 05] 



dtTk = \ Str 



dtR, 



(rf^ 



Ri 



(8) 



where t = ln(fc/A). The function Rk{q^) denotes the 
momentum cutoff. It obeys Rk{q^) — > for k'^/q'^ -^ 0, 
Rkiq^) > for q^/k^ -^ and Rkiq^) -> 00 for fc ^ 
A — )• 00, where fc = A stands for the initial scale in the 
UV. The stability and convergence of the RG flow (Is]) 
is controlled through adapted, optimized choices of the 
momentum cutoff [TT, 39, 49 . Furthermore, F^ denotes 
the second functional derivative of F^; with respect to the 
flelds according to 



(rf) 



ab 






^1_ 



(9) 



where the indices a, b summarize fleld components, inter- 
nal and Lorentz indices as well as coordinates. Note that 
^I^ is merely a collection of fields and not a superfield 



Following the construction in [?7MT| . it is essential 
that the regulator term ASk preserves both the 0{N)- 
symmetry and supersymmetry of the classical theory. Be- 
ing quadratic in the fields it should be the superspace 
integral of ^^Rk{'D'D)Sij^^ . Using the anticommutation 
-2(7'')fc/9^ for the supercovariant 



relation {Vk.'Di} — 
derivatives, we have 



' 1 - \ 2n 

-f>v) = i-ar, 



(10) 



such that a supersymmetric and 0(A^)-invariant regula- 
tor term is the superspace integral of 



$,i?fc(P2?)$' = $, 



ri(-n)-r2(-n)— 1$^ 



Expressed in component fields, we find 



F 



(11) 



(12) 



In momentum space idf^ is replaced by p^ and the bosonic 
and fermionic momentum cutoffs R^ and Rj!^ respectively 
are of the form 



Rl = - {ri + r2f) ® In ■ 



(13) 



Note that the requirements of manifest supersymme- 
try imposes a link between the bosonic and fermionic 
momentum cutoffs, leaving two free functions ri = 
ri(p^/fc^) and r^ = r-iij? jk'^^ at our disposal. Such su- 
persymmetric cutoffs have been introduced for the A'' = 1 
model in two and three dimensions in |551 |3T] . 

There exist no Majorana fermions in three Euclidean 
spacetime dimensions. With respect to the supersym- 
metric 0{N) model we could thus analytically continue 
the fiow equation in Minkowski spacetime to imaginary 
time or alternatively just ignore the fact that the Majo- 
rana condition is not compatible with Lorentz invariance 
in Euclidean spacetime [13]. Both approaches lead to 
identical flow equations in Euclidean spacetime, cf. [ST\ . 



B. Local potential approximation 

Next we turn to the supersymmetric RG flow in the 
local potential approximation. Here, one keeps the lead- 
ing order term in a superderivative expansion such that 
the effective action (with Lorentzian signature) reads 



rfc[$]= Id^x 



1 - / /? 

— ^'VV<^i + 2NWk — 
2 \N 



d^x {d^^'d^^, - i^P'04>, + F^) (14) 



d^x W, 



2N 



where the prime denotes the derivative with respect to 
g/N. The fiow of the renormalized superpotential Wk ( -^ 
in Euclidean space is obtained by projecting the flow Q 
onto the term linear in the auxiliary field F and perform- 
ing a Wick rotation (see appendix^ for its derivation in 
superspace). The function ri acts as IR regulator but 
not as UV-regulator, in contrast to r2 which serves both 
as IR and UV regulator. Thus we use ri as regulator in 
what follows^. Then we find 



d,Wk = - 



d^p 

(27r)3 



dtr2 



N 



,"2)2p2 



M2 



K 



2{g/N)W;! 



N{1 + r2)V + (M^fc + 2{g/N)Wl!y 



(15) 



Similar to the bosonic 0{N) model the flow receives con- 
tributions from the A^ — 1 Goldstone modes (the first 
term) and from the radial mode (second term). 

Next we specify the function r2(p^/fc^). Following [12]- 
[T4] l3T] we choose the optimized regulator function 



r2{p') = 



\P\ 



-l] 






e{k^-p^). 



(16) 



This choice implies dtr2 to vanish identically for p^ > k^, 
and the inverse propagators 



{l + r2)^p^ + X = 



p^ + X for p^ > P 
k'^ +X for p2 < fc2 



become fiat (momentum independent) in the regime 
where the right-hand side of ( |15[ ) is non-vanishing. In 
the LPA this is a solution to the general optimization 
condition of [T^HTl] and is therefore expected to lead 
to improved convergence and stability of the RG fiow. 
Equally important, the momentum integrals in (15) can 
be performed analytically, leading to 



dkWk = 



e 1 



1- 



wi 



Wl+2{-g/N)Wl! 



n' 



87r2 iV fc2 + (w;^ + 2{g/N)W;!y 



(17) 



With given initial condition Wk=A{g/N) = W{g/N) at 
the UV scale A this flow equation uniquely determines 
the superpotential in the infrared limit A: — >■ 0. For A^ = 1 
it reduces to the three-dimensional Wess-Zumino model 
studied in \ST\. 



In order to write the fiow equation (17) in a scale- 
invariant form it is convenient to define a dimension- 
less field variable p as well as a dimensionless superpo- 
tential w and a dimensionless scalar potential v. The 



^ In preliminary studies we did include the regulator ri and got 
almost identical results. 



canonical mass dimension of the fields and potentials are 
[g] = d— 2, [V] = d and \W] = d — 1 in d spacetime 
dimensions. We therefore introduce the dimensionless 
quantities 



N k 



and w{p) — Stt' 



W{jr 



fc2 



(18) 



Note that we have also rescaled an irrelevant numerical 
factor into the potential and the fields. It is understood 
that w is also a function of the RG scale parameter, 
though this is not spelled out explicitly. Similarly, we 
define the dimensionless bosonic potential v as 



v{p) 



87r2 

N 



wa 



= pw'\p), (19) 



where (It]) and ( 18 ) have been used. Thus, by substituting 
(181 into pTJ) we end up with the following flow equation 



for the dimensionless superpotential. 



dtw — pw' +2w 



(1-i)™' j^{uj' + 2pw") 



l + {w' + 2pw"Y ■ ^ ' 



C. Large- A^ limit 

In the large- iV limit the Goldstone modes fully dom- 
inate the dynamics and the contribution of the radial 
mode becomes a subleading effect. It follows that the 
anomalous dimension of the Goldstone modes vanish, as 
no momentum-dependent two-point function exists that 
contribute to the running of the kinetic term of these 
modes to leading order in A''. This is a particular feature 
of the bosonic 0{N) models [J| and their supersymmet- 
ric extensions^. Consequently, the LPA approximation 
becomes exact for A^ — > oo. 

In this limit, the RG equation for the first derivative 
of the superpotential u(p) = w'{p) becomes 



dtu + dpu[l~ p- u^ f{u^)\ = -u 



(21) 



with f{x) = (3 -f x)/{l + xY . We note that the second- 



order partial differential equation ( 20 ) has turned into a 



first-order one in this limit, which is solved analytically 
with the method of characteristics. The first character- 



istic reads 



with 



= const, and the second one is 



p-1 



F{u) 



F{u) = const. 



2 arctan u 



\ + u^ 



(22) 



(23) 



^ For example, Yukawa-type systems may have large anomalous 
dimensions in the large- Af limit |50| . 




FIG. 1: Supersymmetric fixed point solutions p(u*) for all 
fields p and all superfield potentials u* , color-coded by the free 
parameter c (both axes are rescaled as a; — >■ jjt^ for display 
purposes). Thin lines are included to guide the eye, thick lines 
correspond to distinguished values for c (|c| = Q,cl,cp,cm) 
as defined in main text. 



and F'{u) = /(u^). Altogether, we find 



F(u) 



G{ue') 



(24) 



for all p > 0, where the function G(ue*) is determined by 
the boundary conditions for u{p), imposed at the initial 
UV scale k — A. The validity of the solution (24 1 is 
confirmed by direct insertion into ( 21 ). For completeness. 



we also give the RG equation for the bosonic potential. 
Using (19) and (22 1 we obtain 



dtv + 3v — pv' = {v — pv') 



{p + vf 



(25) 



In passing we note that up to minor modifications eq. ( 21 ) 
holds for general spacetime dimensions away from d — i. 
The canonical mass dimension of u is one for all dimen- 
sions and the dependence on spacetime dimensionality, 
therefore, only enters via the field variable leading to the 
replacement of (— p) by (2 — d)p in (21). This modi- 
fies the second characteristic equation whose solution is 
expressed in terms of the hypergcometric function for ar- 
bitrary dimension d ^ 1. Below, we restrict ourselves to 
the case d = 3. 
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FIG. 2: Supersymmetric fixed point solutions it* (p) according to ( 26 1 , covering the entire parameter range for c. With decreasing 
c, fixed point curves rotate counter-clockwise around (p, u*) = (iTOj starting with c = oo where u* = (horizontal line), passing 
through c = (red, dashed-dotted line), completing a rotation of 180° at c = — oo (horizontal line). Further special lines refer 
to \c\ — Cm (blue dashed), \c\ — cp (green, long dashed), \c\ = cl (black, full lines), see main text. Left panel: fixed point 
solutions for all fields (both axes are rescaled as a; — > i£ur\ for display purposes). Right panel: fixed point solutions for physical 
fields in the vicinity of p = 1. 



IV. FIXED POINTS 



A. Supersymmetric fixed points 



Fixed points are the scale-independent solutions of 



through 



(21 1, i. e. solutions satisfying dtu = 0. Besides the Gaus- 
sian fixed-po int solution u* = 0, non-trivial fixed points 
follow from (24 1 in the limit where G(ue*) becomes a 



t-independent constant. The classification of solutions of 



p = 1 + H{u^) + cu^, H{u^,) — u^ F{u^), (26) 



where F[u^) is given by (23 1, then depends only on the 
real parameter c. With |wj]~l [0, oo) and for a fixed c ( 26 ) 
identifies the range of achievable field values. Candidates 
for physical fixed points it* (p) are those solutions which 
extend over all fields p G [0, oo). Fig.Il|and Fig. [2] display 
the entire set of solutions to ( 26 ) for all c. Note that Fig. fTl 
shows the function p(u*) whereas the relation u^{p) is 
displayed in Fig. [2J 

The space of solutions enjoys some internal symme- 
try. Since H{u),) is an even function, solutions only de- 
pend on the absolute value of c, i. e. any solution u*(p) 
with parameter c is equivalent to the reflected solution 
— M*(p) with parameter — c. Both solutions lead to iden- 
tical scalar potentials w* and therefore we may restrict 
our discussion to c > 0. 



(p,w*) = (l,0) 



(27) 



which follows immediately from ( [26| ) due to -ff(O) = 0. 
As can be seen from Fig. [2] the fixed point solutions fall 
into two distinct classes, and solutions in the same class 
show the same global behavior. Depending on the value 
of c the solution u* is either defined for all real p or 
it has a turning point at \pg\ < oo and is only defined 
for p e [ps,oo). In the latter case the solution has two 
branches bifurcating at p = ps- The value of ps will be 
determined below. 

Next we discuss some limiting cases of interest. For 
small u* we conclude from ( 26 1 that 



p-l^cu,+iul + 0{ut). 



(28) 



Hence, the potential is analytical in p — 1 in the vicinity 
of p = 1 for all c, except for c — wher e it becomes 
non-analytical with u* oc ^ p — 1. Eq. ( 28 ) implies that 
all fixed-point solutions have one simple zero at p = 1 
with finite «!(!) except for c — Q where u«(l) diverges. 
Consequently, the scalar fixed-point potentials w* — pu^ 
possess two minima at 



and p = 1 , 



(29) 



We now discuss (26) in more detail. All curves pass 



the first one being a simple zero. The second minimum 
is a double zero for c 7^ and a simple zero for c = 0. 
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FIG. 3: Left panel: Fixed point solutions u« and fixed point potentials «* — pu^ at \c\ = cl showing the two branches w<,u< 
(full lines) and «>,u> (dashed lines). Right panel: The scalar fixed point potential u«(p) as a function of the parameter c 
with Cl (black line), cp (green, long dashed), cm (blue, short dashed) and c — a" cl, a = 2^''^ with n — 1.0,2.3,3.6,4.9 (blue, 
dotted). For cl and cp just one branch is plotted. 



In the largc-u* limit of ( 26 ) , we find 

Thus the asymptotic behavior of m* is given by 
subleading (u* > 0) , 
subleading (u* < 0) . 



C + TT 
P 



(30) 



(31) 



If |c| > Cp, with 



Cp = TT, 



the expansions extend towards p — >■ ±oo, respectively. 
Together with the boundedness of H{u^) we conclude 
that M*(p) is defined for all real p. The expansions corre- 
spond to asymptotically large fields p ^ 1 in the physical 
regime. At \c\ = cp the leading term in (30) vanishes 
and, depending on the sign of c, one of the asymptotic 
solutions is replaced by w* ~ p^^'"^ thus corresponding 
to a small field regime p ^ 1. For \c\ < cp both ex- 
pansions extend towards p — ^ -foo. We conclude that u* 
has, simultaneously, two asymptotic expansions for large 
positive p. This implies that w* displays a loop consisting 
of two branches «< and ti> which coincide at p = cx) and 
possibly at some p — ps < oo where m* has infinite slope. 
The latter condition determines the turning point p^ as 
the simultaneous solution of 



1 



(1 



,2^2 



(33) 



1 



together with ( 26 ) , leading to 

(l+«2)2 

where Ug = u^ {ps)- The degenerate solutions extend over 
the whole physical regime p > 0, provided that ps < 0. 



H{u,) 



(34) 



From ( 33 1 it follows that the equal sign holds for uj = 1 



leading with (34) to \c\ = cl, where 
1 



CL = ;^(7r + 3)« 3.071 



(35) 



For 



c^,, both u< and w> have infinite slope at van- 



(32) ishing field with the non-analytical behavior 



subleading 



(36) 



dur 1 

dp ^~p 

and Us — =f1 for c ~ ±cl (see Fig. |3J left panel). In 
contrast, for cl < \c\ < cp, the behavior at vanishing 
field is analytic. The turning point (ISSJ) exists for small 



< |c| < Cm as long as d'^pldu1\o does not vanish, which 



happens at uf = 3 leading with (34) to \c\ — cm, where 



Cm = 2 




3.177. 



(37) 



We note that cp < cm and conclude that the fixed point 
solutions in the parameter regime cp < \c\ < cm are 
single-valued in the physical regime but multi-valued in 
the non-physical regime p < 0. For all cm < \c\, fixed 
point solutions are single-valued on the entire real axis. 
In Fig. [3] right panel the scalar fixed point potential v^, 
for different values of c is displayed. 



B. Exactly marginal coupling 



C. Line of fixed points 



Next we discuss the physical meaning of the parameter 
c. To this end we employ the polynomial expansion of the 
RG-time dependent superpotential u{t, p) which satisfies 
the flow equation (2f ). For a typical initial condition 
MA = Ti{p — po) there always exists a node po{t) around 
which we can perform a Taylor expansion: 



u(t,p) 



EV«(*)(^-''oW)" 



(38) 



Inserting this ansatz into the flow equation (21) we read 



off the flow equations for po and the couplings r„ entering 
the Taylor expansion 



dtpQ 


= 1- 


dtn 


= 


dtT2 


- ^rl 



Po 



T2 



(39) 
(40) 
(41) 



and similarly to higher order. Several comments are in 
order at this point. Firstly, the running of the vev po{t) 
is independent of all the other local couplings. This prop- 
erty is typical for a supersymmetric flow and has previ- 
ously been observed in [221 EI] ■ The fixed point is ob- 
tained for Po = 1. Secondly, the system of algebraic 
equations describing the i-independent fixed point cou- 
plings can be solved recursively. This leads to fixed point 
couplings r„(ri) for all n > 2 as functions of ti. Insert- 
ing ( 38 ) into the expansion of the scalar field potential 
V — pv? = X]n=2 '^"/"'' (P ~ Po)" and evaluating it on 
the fixed point leads to the fixed point values 



2 


= ^rt 


(42) 


3 


^ 6rf(l-6T2) 


(43) 


3 


= -24rj*(l - 45t2) 


(44) 



In summary, the following picture has emerged. Fixed 
point solutions are characterized by the dimensionless lin- 
ear superfield coupling ti = 1/c in the vicinity of the 
node Po 7^ 0. In the weakly coupled regime 



CP < \c\ 



(46) 



a unique fixed point solution exists covering the whole 
physical domain p > 0. This includes the Gaussian fixed 
point Ti = 0. In the intermediate coupling regime 



Cl < \c\ < cp 



(47) 



two separate fixed point solutions u< and m> exist. The 
former solution has a node at po = 1 whereas the other 
solution has no node, see Fig. |3J left panel. Therefore, 
the corresponding scalar field potentials v^ (?;>) have two 
minima at (29) (one minimum at p = 0). Both are ana- 
lytical functions of p in the vicinity of their global min- 
ima. For \c\ = c/,, the potential becomes non-analytical 
for either of them at p = in a manner reminiscent of the 
Bardeen-Moshe-Bander phenomenon in the purely scalar 
theory [7 . In the strong coupling regime 



\C\ < CL 



(48) 



the theory becomes so strongly coupled that du/dp\p_, 
diverges in the physical regime, and hence no fixed point 
solution exists which extends over all fields. Therefore, 
the supersymmetric 0{N) model displays a line of fixed 
points which bifurcates at \c\ = cp into two fixed points, 
and then terminates at \c\ = cl- 

Finally we note that the solution with c = is closely 
linked to the Wilson-Fisher fixed point in the purely 
bosonic model [33^31] • The precise relation is discussed 
in Sec. |VC] below. 



and similarly to higher order. Clearly, the weak (strong) 
coupling regimes correspond to small (large) A2 and 
hence small (large) ti respectively. Also, on the level 
of the scalar field potential the critical behavior is inde- 
pendent of the sign of ri. Finally, and most importantly, 
the coupling ti remains un-renormalized under the su- 
persymmetric RG flow ( 40 ) . Therefore ri corresponds 



to an exactly marginal coupling, and fixed points can be 
classified according to the value of the linear (dimension- 
less) superfield interaction ti which relates to the free 
parameter c in the analytical solution ( 26 ) as 



1 



V. UNIVERSALITY 
A. Critical exponents 

Fixed point solutions are characterized by universal 
critical scaling exponents. The exponents can be de- 
duced from the RG equations in several ways. Within 
a polynomial approximation up to order n, we expand 
^(p) = J27=i^i(P ~ ^o)V*' in terms of the n -|- 1 cou- 
plings bi{t). From their beta- functions /3i = dtbi, the 
universal exponents follow as the negative of the eigen- 
values 6' of their stability matrix B/ — d(ii/dhj\b=b* 
(45) as Bv^ — —O^v^ with eigenvectors v^ . Using the fiow 



This relation can be shown by inserting expansion ( 38 ) 



into the fixed point equation (26). The presence of the 



exactly marginal coupling ti explains the existence of a 
line of fixed points. 



equation, we find 



= -1,0,1,2,3,. 



(49) 



both numerically and analytically. In fact, the LPA ap- 
proximation has become exact in the large- A^ limit, and 



hence the correct scaHng exponents are achieved to ev- 
ery order in the polynomial approximation. We note that 
this analysis relies on local information of the RG flow 
in the vicinity of m — 0, showing that the scaling (49) 
is achieved mathematically for all < \c\ < oo. Physi- 
cally, however, the analysis is not sensitive to the global 
behavior of the solution, and consequently cannot detect 
that |c| = Ch denotes a physical endpoint. Also, the case 
c = requires special care as an analytical expansion 
about w = is no longer applicable. 



B. Eigenperturbations 

Interestingly, the critical exponents and eigenperturba- 
tions can also be calculated analytically without resort- 
ing to a polynomial expansion. To that end, we consider 
small fluctuations bu about the fixed-point superpoten- 
tial such that u{t,p) = u^,{p) + 5u(t,p). Linearizing the 
flow equation in 6u leads to the fluctuation equation 



dfSu = 



a,-^^^)5u, (50) 



where primes denote a derivative with respect to the 
function's argument. Since the right-hand side is inde- 
pendent of i, the differential equation (50) can be factor- 
ized via separation of variables 5u{t,p) — f{t)g{p) with 



(In/)' 

(Ing)' 



'(lnw*)' + (lnw*0' 



(51) 



where 9 denotes the eigenvalue. Integration leads to the 
exact solution for the linear perturbation of the fixed 
point superpotential 



Su^Ce'^'ul 



(52) 



The allowed range of values for the exponents 9 is deter- 
mined using regularity conditions for the eigenperturba- 
tions. To that end, we recall that the fixed point potential 



w* grows linearly with the field for large p, see (30), and 
hence 6u oc e^*p^^^- Furthermore, 



in the vicinity of the 
node we have (28), which for c ^ leads to a finite u' 



(meaning < i?< oo). We thus find 



du ex e*'* {p - 1) 



9+1 



(53) 



Regularity of the perturbations requires non-negative in- 
teger values for the exponent 9 + 1, reproducing (49). 

Note that this line of reasoning assumes analyticity of 
the perturbation at the node which holds for all c ^ 0. 
For c — 0, u■^, is non-analytical at (27) but ul instead is 



analytical and has a simple zero with finite (u^)'|„,=o- 
Therefore we use ( 52 ) to relate the (regular) fluctuations 
of u^ to ul , leading to 



Su'^Ce'Hulf^^'^'Uuly. 



(54) 



Again, analyticity implies that the exponent {9 + l)/2 is 
a non-negative integer 



-1, 1,3,5,7,. 



(55) 



Here we recognize the universal critical exponents of the 
3d spherical model ^49]. We stress, however, that this 
solution is not a proper fixed point solution in the usual 
sense because it is limited to field values with p > 1. 

Finally we extend the analysis of linear perturbations 
to those of the function u^ and the scalar potential v = 
pv?. We begin with M^ = ul+Sv?' . An analytical solution 
is found by using the identity dv?' = 2u^, du together with 
(52), leading to 



Su^ = 2C e'' ut+^ uL 



(56) 



Note that the degree in u* has increased by one unit. 
Employing the same reasoning as above for c 7^ 0, we 
conclude that the set of available eigenvalues is 



= -2,-1,0,1,2,3,. 



(57) 



Physically, the appearance of the eigenmode with eigen- 
value —2 is due to the mass term squared, a term which 
on dimensional grounds is available in m^ but not in u. 

Finally, using (25), (56) and (21 1, the hnear eigen- 
perturbations about the scalar potential v(t,p) — v^, + 
Sv{t,p) are found as 



5v^2Ce''ul+^u. + u:[l-ulfiul)]} 



(58) 



Close to u* = 0, the term in square brackets reduces to 
1, and the curly bracket becomes u'^ which is finite at 
M* = 0. Therefore regularity of eigenperturbations again 
implies ( |57[ ). In a non-supersymmetric scalar theory the 
potential is not constrained to be of the product form 
(19) and an additional eigenvalue —3 becomes available 



related to redundant shifts of the potential. 

We conclude that supersymmetry is responsible for the 
absence of the redundant eigenvalue —3 in the scalar 
potential, and for relating its two relevant eigendirec- 
tions with eigenvalues —1 and —2 with the sole relevant 
eigendirection with eigenvalue —1 of the derivative of the 
superpotential. 



C. Wilson-Fisher fixed point 

It is interesting to clarify how the supersymmetric 
model and its fixed points fall back onto those of the 3d 
non-supersymmetric scalar theory in the same approxi- 
mation [37K39] . To that end, we consider the 4d super- 
symmetric 0{N) at finite temperature. The tempera- 
ture is implemented using the imaginary time formalism 
which on the level of the flow equation amounts to the 
replacements [101 US [15] 



dqo 
27r 



fiqo) ^ T 



00 

E 

71 — — oc 



f{qo = 27rc„T) , 



(59) 



Here 27rc„T denotes the nth Matsubara frequency with 
Cn = n for bosons and c„ = n + 2 ^r fermions. The tem- 
perature imposes periodic (anti-periodic) boundary con- 
ditions for bosons (fermions) and, consequently, softly 
breaks global supersymmetry. Within a derivative ex- 
pansion the relevant momentum integrals are performed 
analytically using the four-dimensional version of ( 15 ) 



(161. 



together with ( 59 ) and the optimized momentum cutoff 



We are interested in the large-scale behavior k/T -^ 0. 
Due to ( 59 ) , all fermions and bosons with a non- vanishing 



Matsubara mass will decouple from the system, except 
for the bosonic zero mode. In this limit, the Ad super- 
symmetric model undergoes a dimensional reduction to 
a 3d non-supersymmetric theory where all fermions have 
decoupled. In the large- A^ limit, the RG equation for the 
potential of the remaining bosonic zero mode in LPA is 
given by 




FIG. 4: The Wilson-Fisher fixed point solution Zt.{p) of (60 1 



dfZ 



-2z + pz' 



1 



(i + zy 



(60) 



where z is related to the scalar field potential by v{p) = 
pz{p). The key difference to the supersymmetric sys- 
tem studied previously is that the function z is no longer 
constrained to be the square of a superpotential deriva- 
tive w' . Relaxing this constraint allows for an additional 
fixed point solution, which follows from integrating ( 60 ) 
analytically. The general solution reads 



1 



V^ 



l + z 



2 arctan \/z = B{z e 



2t\ 



(61) 



where B{ze'^*) is fixed through initial conditions. The 
solution for negative z is found by analytical continua- 
tion. In particular, (60) has a Wilson-Fisher fixed point 
solution z* 7^ with z(/9 = 1) = corresponding to (61) 
with B = 0. The solution extends over all p with one 
unstable direction, see Fig. l4J The eigenperturbations 
Sz are found analytically leading to (54) with 
6z and ui — >■ 



z ^ z, 

the replacements Su 

universal eigenvalues are identical and given by ( 55 ) 



Hence, the 



The similarities and differences between the Wilson- 
Fisher fixed point solution of the purely scalar theory 
and the c — 'would-be' Wilson-Fisher fixed point of 
the supersymmetric partner theory can also be appre- 
ciated from the behavior at small and large fields. In 
fact, for p > 1, z.^,{p) is positive and related to the real 
superpotential by 



z*{p) ^w'^ipf 



In turn, z^,[p) is negative for all p < 1. Interestingly, 
this solution is still visible in the supersymmetric theory 
where it corresponds to a purely imaginary "superpoten- 
tial" with 



Hence, provided that a purely imaginary superpotential 
is meaningful in the supersymmetric theory, the c = so- 
lution can be extended to a valid supersymmetric Wilson- 
Fisher fixed point for all p. However, the structure of 
the Lagrangian imposed by supersymmetry implies that 
the field-dependent fermion mass term is proportional 
to w'^ and the Yukawa-type fermion-boson interaction 
proportional to w'l all become purely imaginary. Most 
importantly, a purely imaginary w'^ for small fields im- 
plies that the scalar potential obeys v*(p) = pw'^ < 
for all fields within < p < 1. Unbroken global super- 
symmetry requires that the dimensionful Vk{g) remains 
positive for all fields and scales. In the infrared limit 
fc — >■ 0, re-inserting powers of fc, the dimensionful po- 
tential approaches V{g) = GAir'^g^/N'^ > 0. Hence, our 
results state that this potential can be approached arbi- 
trarily close from within a phase with 0{N) symmetry 
and global supersymmetry. 



VI. CONCLUSIONS 



We have studied fixed points of supersymmetric 0{N) 
symmetric Wess-Zumino models in the limit of many 
components TV —> 00 in three dimensions with the help 
of the renormalization group. We have solved the the- 
ory analytically, showing that it displays a line of non- 
trivial fixed points solely parametrized by the exactly 
marginal linear superfield coupling. The fixed points are 
non-Gaussian, yet they display Gaussian exponents sim- 
ilar to the line of fixed points observed in the bosonic 
(62) (0^)"^ theory. The line of fixed points contains the Gaus- 
sian fixed point and therefore all fixed-points are continu- 
ously linked to the Gaussian one. With increasing super- 
field coupling, the line of fixed points bifurcates into two 
fixed point solutions, both of which terminate at a crit- 
ical coupling ( 35 ) below which no fixed point solutions 



w'Ap) ^±i\/~z*{p)' 



(63) 



exist which extends over all physical fields. One of these 
solutions has its minimum at po = the other at po 7^ 0. 
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Interestingly, remnants of the non-Gaussian scaling ex- 



ponents of the 3d spherical model ( 55 ) become visible for 
asymptotically large superfield coupling. However, the 
fixed point solution does not extend over all fields in the 
supersymmetric case, except if the superfield potential 
becomes purely imaginary for small fields. 

From a structural point of view, the main impact of 
global supersymmetry on the critical behavior in com- 
parison with the purely scalar theory is summarized 
as follows. Firstly, for unbroken global supersymmetry 
the scalar potential has its minimum at vanishing field. 
Hence the irrelevant eigenmode with eigenvalue —3 corre- 
sponding to overall shifts in the potential is absent from 
the supersymmetric eigenvalue spectrum. Secondly, the 
quartic and sextic coupling of the scalar potential are 
no longer independent. Hence, in the supersymmetric 
theory criticality is achieved by tuning only one param- 
eter as opposed to the tuning of two parameters in the 
corresponding purely bosonic theory. This is refiected 
in the sole negative eigenvalue for u as opposed to the 
two negative eigenvalues for both u^ and v. Finally, at 
the coupling \c\ = cl (35) the supersymmetric model 



shares similarities with the Bardeen-Moshe-Bander phe- 
nomenon in the bosonic theory [8j. The logarithmic sin- 
gularity observed in J8] is superseded by a square-root 
behavior in the supersymmetric case, a difference which 
can be traced back to the underlying regularizations. 

The fixed point solutions discussed in this paper de- 
scribe the phase transition for the breaking of the 0{N) 
symmetry. Analyzing the pattern of symmetry breaking 
and the phase transition between symmetric and bro- 
ken phases in more detail, and relating our findings with 
earlier studies based on gap equations is deferred to an 
upcoming publication. Furthermore, stepping back to 
finite N we expect modifications to the above picture, 
both within the local potential approximation studied 
here and to higher order in the derivative expansion. For 
example, it is known that the A^ = 1 model displays a 
superscaling relation linking the unstable direction with 



the anomalous dimension [551 ISI] : a behaviour which is 
quite different from the Ising universality class [24]. It 
will thus be interesting to see how these patterns gener- 
alize for supersymmetric 0{N) models with generic N . 
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Appendix A: Conventions 

Relevant symmetry relations and Fierz identities for 
Majorana spinors are ^x — X^j ^"f^x — — X7''^ and 
QkQi = —■^{()0)lki- One of the main features of the action 
is its invariance under supersymmetry transformations. 
The latter are characterized by the supersymmetry vari- 
ations 5^^^, generated by the J\f — I fermionic generator 
Q. We have 

S,<i>\x) = iekQk<^\x) with Qk - -tde^ ~lkAd^.. (Al) 

Thus, (JATJ) leads to the supersymmetry variations 

S(l)' = e^\ (5?/'' = {F' + i$(f)')e and SF' = ie$^P'' (A2) 

of the component fields. The anticommuting sector of 
the superalgebra is given by the anticommutator of two 
supercharges 



{Q,,Qi}^H,d^. 



(A3) 



Appendix B: Superspace 

Following [57] we consider the action of the three-dimensional supersymmetric 0{N) model in the local potential 
approximation 



Vkm = jd^x 



^(-l$^A^$,: + 2iy,(i?)), 



(Bl) 



where R — ^(E'^^i, K — ^{VD — VD) and i — 1, ..., A^. We derive the fiow equation in the superspace M"^'^ with 
coordinates z — (x, 01,62)- Furthermore, we introduce the abbreviation J dz = J d^xd9id92/{2i). In Minkowski 
spacetime [i31_, the Wetterich equation in superspace may be written in the form 



dtTk = ^ ldzdz'{dtRk)mu{z,z'){Gk)nm{z\z), t^H^/K^), 



(B2) 



where {Rk)mn represents a supersymmetric regulator term and {Gk)nm the connected Green's function. According 
to [571 EI] J we now choose a general regulator term quadratic in the superfields $* and diagonal with respect to the 
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field indices: 



Notice that this regulator conserves both the 0{N) symmetry and supersymmetry. The functional derivative with 

f^ / dz $^(z) = j;. with g;{|)= 2*5;. 



respect to a superfield is chosen according to the conventions ^^f.-^ J dz $*(z) = (5* with ^^j-.lj = 2i <5* 5{x — x) (5(02 " 



62) S{9i — 61) = SjS{z — z). Thus, the second functional derivative of the effective average action with respect to the 
superfields reads 



ri'L(^. ^') ^ Ach^^fe— ^ = [{-K + 2W;,{R))6,,^ + 2<(i?)cl>„ci>„] {z)S{z - z'). (B4) 



5$"(z)^'''(5$"(z') 

Similarly, the second functional derivative AS), (z, z') of the regulator term is given by 

(i?fe)„„,(z, z') - [2ri - r2K] (z) (5„„ 5{z ~ z'). (B5) 

Now we assume the superfields to be constant, i.e. 9^<i>*(x, 0) = 0, such that the regulator functions as well as the 
wave operator may be simply written in momentum space. However, note that the wave operator K still contains 
derivatives with respect to the Grassmann coordinates and thus acts on the adjacent delta functions. Hence, the flow 
of the effective average action may be written as 

dtVk - \ I dzdz' {dtRk)mn{z,z'){Tf +Rk);,l{z',z) 



I 



dOi d92 d0[ de'2 f d^p 



- I d^x'^^^^'^^^^J ^^i2dtn-dtr2K)ip,9,,d2)S,nn2tSie2-9'^)6iei-9[)x 

[i~hK^p^e[,e',) + 2W')5„™ + 2W"$„$™] ''2i 5[9'^ - ^2) 5{9[ - 9^). (B6) 

We have thereby introduced the notation W^R) = Wl{R) + ri, h= 1 + r2. The inverse of the N x iV-matrix 

{M)-l = {~hK + 2W')'5„™ + 2W"$„$„, (B7) 

is given by 

, . ^ j-hK + 2W')Sr,m + 2W"i^^S^^ - $„$„) 

^ ^"" 4{h?p^+W'(W' + 2W"R)-hK{W' + W"R))' ^ ' 



where we have used the relation K'^{p) = 4p-^ resulting from the action of K{p) = —dedg — {de^9) — {9^dg) — p 



on an arbitrary superfield. In order to eliminate the wave operator K in the denominator of ( B8 ) , we multiply both 
the numerator and the denominator with [{h'^p'^ + >V'(W + 2W"R}) + hK{W' + W"R)] and use again K'^{p) = 4p2. 
Thus we get 

^ '""" 4(/iV- W2)(/iV-(>V' + 2W"i?)2) 

.Snrnih^P^ " W'^) - 2W" {W + W"i?)($2^„„, - $„$„)) 



-hK- 

4(/i^p^ - W'^)[h^p^ - {W + 2W"R)^) 

(B9) 



4(/iV - W^){h^p^ - (W + 2W"E)2) 

n 



with n = A{h'^p^ - W2)(/iV _ (yy/ _^ 2W"i?)2). For 

iGk)nm(j>,e[^e,,e',^e,) = (M)„™(p, 0^, 0^) 2t S{9'2 - 92)S{9[ - 9,) (BIO) 

to be the Green's function it has to fulfill the defining relation 

dz iGk)rnn{~Z, z){r'-k^ + Rk)np{z, z') = 6{~Z - z')(5™p. (Bll) 



This can be shown by directly inserting the explicit expressions on the left hand side and working out the contributions 
to different orders in K. 
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The flow equation is calculated by inserting the regulator ( |B5[ ) as well as the propagator (BIO I into eq. (B6). Note 
that the regulator {Rk)mn oc 5mn is diagonal with respect to the field indices. Hence, we simply have evaluate the 
trace over the Green's function {Gk)mm- This yields 






d^x 



d'^p 

(2^ 



{2dtn - dtr2K){p, 9^,02) 6(92 - 0^) <5(0i - e[) 



h^p^{NW' + 2W"R) - W'(W + 2W"R){NW + 2{N - l)W"R) 



^hK 



4:{h^p^ - W'^){h^p^ - (W + 2>V"i?)2) 
JV(/iV - W) - m - 1)W"{W' + W"R)R \ _ , _ 

4(/iV-W'2)(/iV-(W' + 2W"i?)2) J^^'^^ "^'^^^^ "^>- 



(B12) 



Now, only terms linear in K contribute to the flow of Ffc after having integrated out the Grassmann variables. Those 
contributing terms lead to a multiplying factor of 2i. Thus the flow equation (B12| simplifies to 



dtTk 



d^xd9idd2dtWk{R) 



- Id^xdeide2 



d^ A^ _ jdtnh - dtr2W') ^ dtrih-dtr2iW' + 2W"R) 



(2^) 



/j2p2 _ ^;,2 



/j2p2 _ (yy. ^ 2W"i?)2 



(B13) 



Performing a Wick rotation of the zeroth component of the momentum, i.e. p" — >■ ip%, p^ — > — p|;, we obtain the 
Euclidean version of the flow equation (B13). Thus, the resulting flow equation in superspace reads 



dtWkiR) = - 

X,tli,62 ^ •'XM1M2 



(PpE 

(27r)3 



(iV-1) 



{dtnh - dtr2W') dtnh - dtr2{W' + 2W"R) 



h^pl + W2 



h^pl 



(W + 2W"i?)2 



(B14) 



Notice, that the truncation p4| involved a superpotential of the form 2NWk{R/N) instead of 2Wk.{R)- The cor- 
responding flow equation may be easily derived from the above result by performing the substitution Wk{R) — > 
NWkiR/N) in (|B14|). This yields the final result 



/ 



dtWkiR/N) = 



d^PE f iN-l){dtnh-dtr2W') 1 dtnh-dMW + 2W"R/N) \ 

3 \ N /.2„2_ _^^,2 + JSf h2p2^ + (y^f _^ 2'W"R/NY J ' ^ ' 



(2^) 



h^pl 



r 



[1] Kenneth G. Wilson. Renormalization group and critical 
phenomena. 1. Renormalization group and the Kadanoff 
scaling picture. Phys. Rev., B4:3174-3183, 1971. 

[2] Kenneth G. Wilson. Renormalization group and critical 
phenomena. 2. Phase space cell analysis of critical behav- 
ior. Phys. Rev., B4:3184-3205, 1971. 

[3] Jean Zinn-Justin. Quantum Field Theory and Critical 
Phenomena. Oxford University Press, third edition edi- 
tion, 1996. 

[4] S. Weinberg. Ultraviolet divergences in quantum theories 
of gravitation, in: General Relativity: An Einstein cen- 
tenary survey, Eds. S. W. Hawking and W. Israel. Cam- 
bridge University Press, first edition edition, 1979. 790- 
831. 

[5] Daniel F. Litim. Renormalisation group and the Planck 
scale. Phil. Trans. R. Soc, A 369:2759-2778, 2011. 
arXiv:1102.4624 [hep-th]. 

[6] Andrea Pelissetto and Ettore Vicari. Critical phenomena 
and renormalization group theory. Phys.Rept., 368:549- 
727, 2002. cond-mat/0012164 

[7] William A. Bardeen, Moshe Moshe, and Myron Bander. 
Spontaneous Breaking of Scale Invariance and the Ultra- 



violet Fixed Point in 0(n) Symmetric (phi^ in Three- 
Dimensions) Theory Phys. Rev. Lett, 52:1188, 1984. 
[8] Francois David, David A. Kessler, and Herbert Neu- 
berger. The Bardeen-Moshe-Bander Fixed Point and 
the Ultraviolet Triviality of phi**6 IN three-dimensions. 
Phys. Rev. Lett, 53:2071, 1984. 
[9] Joseph Polchinski. Renormalization and Effective La- 
grangians. Nucl. Phys., B231:269-295, 1984. 

[10] Jurgen Berges, Nikolaos Tetradis, and Christof Wet- 
terich. Non-perturbative renormalization flow in quan- 
tum field theory and statistical physics. Phys. Rept., 
363:223-386, 2002. hep-ph/0005122 

[11] C. Bagnuls and C. Bervillier. Exact renormalization 
group equations: An introductory review. Phys. Rept., 
348:91, 



, 2001. |hep-th / 0002034, 
[12] Daniel F. Litim. Optimised renormalisation group flows. 

Phys. Rev., D64:105007, 2001. hep-th/0103195 
[13] Daniel F. Litim. Optimisation of the exact renormal- 
isation group. Phys. Lett, B486:92-99, 2000. hep-. 



Jh/0005245 
TI41 Daniel F. Litim. Mind the gap. 
A16:2081-2088, 2001. hep-th/0 104221 



Int.J.Mod.Phys., 



13 



t 



[15] Daniel F. Litim. Wilsonian flow equation and thermal [33 
field theory. 1998. ^hep-ph/9811272 

[16] Daniel F. Litim and Jan M. Pawlowski. On gauge invari- 
ant Wilsonian flows. 1998. hep-th/9901063 [34 

[17] K. Aoki. Introduction to the nonperturbative renormal- 
ization group and its recent applications. Int. J. Mod. 
Phys., B14:1249-1326, 2000. [35 

[18] B. Delamotte, D. Mouhanna, and M. Tissier. Non- 
perturbative renormalization group approach to frus- 

trated magnets. Phys. Rev., B69:134413, 2004. | cond- | 

mat/0309101, [36 

Jan M. Pawlowski. Aspects of the functional renormal 
jsation group. Annals Phys., 322:2831-2915, 2007. hep 
th/0512261, [37 

Holger Gies. Introduction to the functional RG and ap- 
plications to gauge theories. 2006. hep-ph/0611146 
Martin Renter and Frank Saueressig. Functional Renor- [38 
malization Group Equations, Asymptotic Safety, and 
Quantum Einstein Gravity. 2007. arXiv:0708.1317 [hep- [39 
th]. 

Hidenori Sonoda. The Exact Renormalization Group - 
renormalization theory revisited -. 2007. |arXiv:0710. 1662| [40 
[hep-th]. 

Oliver J. Rosten. Fundamentals of the Exact Renormal- 
ization Group. 2010. arXiv:1003. 1366 [hep-th]. 
Daniel F. Litim and Dario Zappala. Ising exponents 



t 



20] 

[21] 



[22] 

[23] 
[24] 

[25; 

[26; 

[27 

[28; 
[29] 

[30] 

[3i; 

[32; 



[41 

from the functional renormalisation group. Phys. Rev., 
D83:085009, 2011. arXiv:1009.1948 [hep-th]. 
F. Vian. Supersymmetric gauge theories in the exact [42 
renormalization group approach. 1998. hep-th/9811055 
M. Bonini and F. Vian. Wilson renormalization group 
for supersymmetric gauge theories and gauge anomalies. [43 
Nud. Phys., B532:473-497, 1998. hep-th/9802196 
Franziska Synatschke, Georg Bergner, Holger Gies, and 
Andreas Wipf. Flow Equation for Supersymmetric Quan- [44 
tum Mechanics. JHEP, 03:028, 2009. arXiv:0809.4396 
[hep-th]. 

Holger Gies, Franziska Synatschke, and Andreas Wipf. 
Supersymmetry breaking as a quantum phase transition. [45 
Phys. Rev., D80:101701, 2009. arXiv:0906.5492 [hep-th]. 
Franziska Synatschke, Holger Gies, and Andreas Wipf. [46 
Phase Diagram and Fixed-Point Structure of two di- 
mensional N=l Wess-Zumino Models. Phys. Rev., 
D80:085007, 2009. arXiv:0907.4229 [hep-th]. [47 

Franziska Synatschke-Czerwonka, Thomas Fischbacher, 
and Georg Bergner. The two dimensional N=(2,2) 
Wess-Zumino Model in the Functional Renormaliza- 
tion Group Approach. Phys. Rev., D82:085003, 2010. [48 
'arXiv:1006.1823 [hep-th]. 

Franziska Synatschke, Jens Braun, and Andreas Wipf. [49 
N=l Wess Zumino Model in d=3 at zero and flnite tem- 
perature. Phys. Rev., D81:125001, 2010. arXiv:1001.2399| [50 
[hep-th]. 

Franziska Synatschke, Holger Gies, and Andreas Wipf. 
The Phase Diagram for Wess-Zumino Models. AIP Conf. 
Proc, 1200:1097-1100, 2010. arXiv:0909^4l89J [hep-th]. 



Sven Falkenberg and Bodo Geyer. Effective average ac- 
tion in N = 1 s uper- Yang-Mills theory. Phys. Rev., 
D58:085004, 1998. |hep-th/9802Tl3| 

Oliver J. Rosten. On the Renormalization of Theo- 
ries^ a Scalar Chiral Superfleld. JHEP, 03:004, 2010. 
arXiv:0808.2150 [hep-th]. 

Hidenori Sonoda and Kayhan Ulker. An elementary 
proof of the non-renormalization theorem for the Wess- 
Zumino model. Prog. Theor. Phys., 123:989-1002, 2010. 
arXiv:0909.2976 [hep-th]. 

Hidenori Sonoda and Kayhan Ulker. Construction of a 
Wilson action for the Wess-Zumino model. Prog. Theor. 
Phys., 120:197-230, 2008. arXiv:0804.1072 [hep-th]. 
N. Tetradis and D. F. Litim. Analytical Solutions of 
Exact Renormalization Group Equations. Nucl. Phys., 
B464:492-511, 1996. hep-th/95 12073 
D. Litim and N. Tetradis. Approximate solutions of exact 
renormalization group equations. 1995. hep-th/9501042 
Daniel F. Litim. Critical exponents from optimized renor- 
mafization group flows. Nucl.Phys., 6631:128-158, 2002. 
hep-th/0203006 

John F. Dawson, Bogdan Mihaila, Per Berglund, and 
Fred Cooper. Supersymmetric approximations to the 3D 
supersymmetric 0(N) model. Phys. Rev., D73:016007, 
2006. hep-ph/0512016 

William A. Bardeen, Kiyoshi Higashijima, and Moshe 
Moshe. Spontaneous Breaking of Scale Invariance in a 
Supersymmetric Model. Nud. Phys., B250:437, 1985. 
Moshe Moshe and Jean Zinn- Justin. Quantum field the- 
ory in the large N limit: A review. Phys. Rept., 385:69- 
228, 2003. hep-th/0306133 

Moshe Moshe and Jean Zinn- Justin. Phase structure 
of supersymmetric models at finite temperature. Nud. 
Phys., B648: 131-160, 2003. hep-th/0209045 
J. Feinberg, M. Moshe, Michael Smolkin, and J. Zinn- 
Justin. Spontaneous breaking of scale invariance and su- 
persymmetric models at finite temperature. Int. J. Mod. 
Phys., A20:4475-4483, 2005. 

Tsuneo Suzuki. Three-Dimensional 0(N) Model with 
Fermi and Scalar Fields. Phys. Rev., D32:1017, 1985. 
Ragnheidur Gudmundsdottir and Gunnar Rydnell. On a 
supersymmetric Version of (phi**2)**3 in three- dimen- 
sions Theory. Nud. Phys., B254:593, 1985. 
Yoshimi Matsubara, Tsuneo Suzuki, Hisashi Yamamoto, 
and Ichiro Yotsuyanagi. On a Phase with Sponta- 
neously Broken Scale Invariance in Three-Dimensional 
0(N) Models. Prog. Theor. Phys., 78:760, 1987. 
Christof Wetterich. Exact evolution equation for the ef- 
fective potential. Phys. Lett, B301:90-94, 1993. 
Daniel F. Litim. Derivative expansion and renormaliza- 
tion group fiows. JHEP, 0111:059, 2001. hep-th/0111159 
Jens Braun, Holger Gies, and Daniel D. Scherer. Asymp- 
totic safety: a simple example. Phys. Rev., D83:085012, 
2011. arXiv: 1011.1456| [hep-th]. 



